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Abstract 
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The x- and Q 2 -dependences of the Gottfried sum rule Sg(x,Q 2 ) based on the 
^ experimental data on proton and deuteron structure functions are studied. The 

dependence of Sg{x, Q 2 ) on Q 2 for low x points to flavour asymmetry of sea quark 
Qh! distributions. For the first time, the coefficients ci, C2 of the expansion of Sg{x, Q 2 ) 

in as/ft up to second order are evaluated through the phenomenological analysis 
of NMC and HI data. It is found that c\ is negative while C2 is positive. The 
obtained result is in disagreement with QCD predictions for these coefficients. We 
suggest measuring the low x-dependence of the F 2 D deuteron structure function at 



2 

2 



HERA in order to study the Q -dependence of the Gottfried sum rule 



1 Introduction 



The experimental data on the proton and neutron structure functions (SF) are of 
great interest for verification of the theory of strong interaction, QCD. The relevant 
information can be used to extract the spin-dependent and spin-independent parton 
distributions, to estimate nonperturbative effects, to verify nucleon models and sum 
rules such as the Gottfried |TJ, Bjorken [Q, Ellis- Jaffe |JJ, Gross-Llewellyn Smith || and 
Adler0 ones. 

New data on the deuteron SF F 2 D obtained at CERN, SLAC and Fermilab @-[|] 
stimulated great interest in theoretical studies of the deuteron structure. 

The deuteron is an excellent neutron target and therefore the neutron structure func- 
tions F£, g™ were usually extracted from the experimentally known proton and deuteron 
structure functions. 

The extraction procedure of the neutron SF from deuteron and proton data is ambigu- 
ous and, therefore, the estimate of nuclear effects in the deuteron is extremely important 
not only to obtain new information on F£ but also to verify deuteron models and to 
perform a common QCD analysis of experimental data. 

The Gottfried sum rule was verified by the NMC Collaboration || and the value of 
Sq = 0.240 ± 0.016 was found to be below the parton model prediction. To study the 
sum rule in detail, the experimental data on the deuteron structure function F 2 D at low 
x and large Q 2 are necessary. Such measurements are possible to perform at HERA JTH . 

In the present paper, the phenomenological analysis of x and Q 2 dependences of 



the S G (x,Q 2 ) Gottfried sum rule is based on the NMC 0, g, HI [|TTJ and ZEUS p 



parametrization of the proton structure function F 2 (x, Q 2 ). It is shown that the available 
experimental data on F 2 and F 2 allow one to estimate the as correction to Sg{Q 2 )- 
The expansion coefficients of Sg(Q 2 ) up to order Q(a 2 ) are estimated. It is found 
that the obtained results are in disagreement with the QCD predictions. To clarify the 
discrepancy, it was proposed to measure the F 2 D (x,Q 2 ) deuteron structure function at 
low x and high Q 2 in order to extract precisely the coefficients of a^-expansion of the 
Gottfried sum rule. 



2 Deep-Inelastic Scattering on Deuteron 

The cross section of deep - inelastic lepton - deuteron scattering in the one - photon 
approximation is expressed via the imaginary part of the forward scattering amplitude 
of the virtual photon or W-boson on the deuteron - W®. The latter is related to 
the deuteron spin-dependent - g\ 2 {y-,Q 2 ) and spin-independent - Ff 2 3 (V, Q 2 ) structure 
functions as follows 

= ~(g^ - W<? 2 ) • + (p„ - qM/q 2 )(p„ - qu(pq)/q 2 ) • F 2 D /u 

+ie^q a {s P 9?lv + [AOP) ~ V P (sq)}M- l g»/u 2 } + ie^ q a p P ■ F 3 D /u. (1) 



Here q,p are momenta of the photon and deuteron; M is the deuteron mass; v = (pq); 
the 4-vector s a describes the deuteron spin. The symmetric part of the deuteron tensor 
W£ can be written as W° = W$ ■ pJ2 



W a(3 
' ' [11/ 
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5(m 2 - k 2 )6(k )6(p + - k+) Sp{w" ■ r(h) -(m + k)- ^(h)}. (2) 



' [IV 



Here the ^-function and light-cone variables - (k±,k±) are used. The tensor p K a J is the 
symmetric part of the deuteron polarization density matrix. The antisymmetric part 
of the deuteron tensor W® v is expressed in a form similar to (2). The procedure to 
construct the relativistic deuteron wave function (RDWF) ip a was proposed and RDWF 
was obtained in |[13 



3 Deuteron and Neutron Structure Functions 

The deuteron SF in the light-cone variables is expressed as follows 

F 2 D (a, Q 2 ) = f 1 dx d 2 k ± p(x, k ± ) ■ F?(a/x, Q 2 ). (3) 

The nucleon SF = (F% + F£)/2 is defined by the proton and neutron ones. The func- 
tion p(x, k±) describes the probability that the active nucleon carries away the fraction 
of the deuteron momentum x = ki + /p + and the transverse momentum k± in the infinite 
momentum frame. It is expressed via the RDWF 



p(x, k ± ) oc Spilth) ■ (m + h) ■ ^(h) ■ q/u ■ pg}. (4) 
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The nuclear effect in the deuteron is described by the ratio R F = Fr P/F% ■ It was 
shown in [1^| that the effect of relativistic Fermi motion grows with x and the ratio R^ N 



reaches 6% at x ~ 0.7. The dependence of the ratio Rp on x resembles the nuclear 
EMC effect and it is practically independent of Q 2 . Using the universal behaviour of the 
ratio it is possible to extract the neutron SF F£ 

F?(x, Q 2 ) = 2- [R^ix)] 1 ■ F 2 D (x, Q 2 ) - F%(x, Q 2 ). (5) 

We would like to note that for a reliable estimate of other contributions to F® from 
nuclear effects such as nuclear shadowing, meson exchanges etc., the data at low x such 
as the E665 data [[| but with smaller experimental errors are required. 

4 Gottfried Sum Rule 

The extracted neutron SF can be used to verify the Gottfried sum rule PJ: 

" [F$(x) - F%(x)]dx/x = 1/3. (6) 
o 



The Gottfried integral as a function of x and Q 2 is defined as follows 



S G (x,Q 2 ) = [\Fi(y,Q 2 )-F?(y,Q 2 )]dy/y. (7) 

J X 

To verify the sum rule, not only the x-dependence of SF's also Q 2 -dependence in 
a wide kinematical range are necessary. The realistic comparison of the experimental 
results with theoretical predictions based on QCD is more argued at high Q 2 . 

As has been reported in ||, the value of Sq at Q 2 = 4.0 (GeV/c) 2 obtained from the 
measurements of F 2 and F 2 is considerably below the value of the naive quark-parton 
model equal to 1/3: Sq = 0.240 ± 0.016. This result in the parton model is usually 
interpreted as the violation of the isospin symmetric sea. As will be shown later, the 
asymmetry violation demonstrates the strong Q 2 -dependence. 

The QCD corrections of order 0(a s ) \TE\ and 0(a 2 ) |Tj| are estimated for the Got- 
tfried sum rule in the case of flavor-symmetric sea u = d: 

Sg(Q 2 ) = \ ■ (1 + c x ■ (a s M + c 2 ■ (asM 2 ). (8) 

The coefficients Ci,c 2 are equal to 0.036, 0.72 for rif = 3 and 0.038, 0.55 for n/ = 4, 
respectively. Thus, the coefficients Ci, c 2 are found to be positive and relatively small, 
and as mentioned in ||16|| , the QCD corrections cannot explain the deviation of the 



theoretical prediction from the experimental result of the NMC collaboration without 
the assumption that the light quark sea is flavor asymmetric. 

5 Procedure to Extract Neutron Structure Func- 
tion 

The method to extract F 2 {x, Q 2 ) from proton and deuteron experimental data was 
proposed and realized in fl4l . 

The procedure includes the items 

• experimental data on the ratio Rp^ p = F 2 D / F 2 and structure functions F 2} F 2 

• the relativistic deuteron model [0 

• the choice of the parametrization of the F 2 neutron structure function and the 
determination of free parameters to describe the ratio Rp^ p = F 2 D / F 2 

• the comparison of the absolute values of the experimental and theoretical structure 
function F 2 D (x,Q 2 ) 

• the determination of the ratio Rp^ N = F 2 / F 2 describing the nuclear effect in the 
deuteron 

• the extraction of the neutron structure function ^from experimental data using the 
formula 

F?(x, Q 2 )=2- [R^Y 1 ■ F 2 D (x, Q 2 ) - F|(x, Q 2 ) 



In fTJ the NMC data g |7J on the ratio R 



D,p = Ff/F%, F% and the relativistic 
deuteron model were used to extract the neutron SF F 2 n . It was shown that the cal- 
culated results for ratio Rp /P and F®(x,Q 



and F^(x,Q 2 ) are in good agreement with the available 
experimental data [§, [7], [H], [I8[] . 

Thus, the conclusion was made that the extraction procedure proposed for F£(x, Q 2 ) 
is self-consistent because it provides a good description of higher statistics experimental 
data on the ratio R^ p and F® over a wide kinematic range of x and Q 2 . 

We would like to emphasize that the nuclear effect of Fermi motion was only used in 
|i~4f . The shadowing effect |L9| should be also included into the procedure if the deuteron 
data in the low x- and Q 2 - range are used ( for example, E665 data). In that case, the 
factor R D I N should be corrected at low x. 



6 Results and Discussion 

Figure l(a,c,e) shows the dependence of the Gottfried integral Sg(x,Q 2 ) on x and 
Q 2 . The parametrization of the proton structure function F%(x, Q 2 ) for NMC and HI 
data are taken from J7|, |8|, O. The NMC parametrization was used in |TJJ for calculating 
the deuteron structure function -F 2 D (x, Q 2 ) and a good agreement with experimental data 
SLAC, BCDMS, NMC was obtained both for a low and high x range. It is assumed [TJJ 
that the ratio f(x) = F£(x, Q 2 ) / F^x, Q 2 ) — > 1 as x — > 0. The parametrization of the 
f{x) function was obtained and will be used later for the determination of the neutron 
SF from the HI and ZEUS proton structure function parametrization. 

One can see from Figure l(a,c,e) that the x- and Q 2 - dependences of Sq(x, Q 2 ) for the 
parametrizations are similar. We would like to note that there is the crossover point Xq 
( in particular, x ~ 0.01 for the HI parametrization) separating two ranges: one - with 
decreasing Sa(x, Q 2 ) and the other - with increasing Sg(x, Q 2 ) with Q 2 , respectively. A 
similar dependence of the Sgls{x, Q 2 ) Gross-Llewellyn Smith integral on Q 2 with the 
crossover point xq =~ 10~ 2 is predicted in [p0~| . Note that the main part of the Sg(x, Q 2 ) 
integral is given by integration over the interval of relatively large x ( 10~ 2 < x < 1.) 
and determines the first term in expansion (8). 

Figure l(b,d,f) shows the dependence of the derivative dSc(x, Q 2 )/das on x and Q 2 . 
The 3-loop perturbative QCD expression of as with rif = 4 is used. The derivative 
grows up to x ~ 0.1 — 0.2 then decreases and changes in sign for Q 2 > 100 GeV 2 for 
NMC and even for smaller Q 2 for HI parametrization. It should be stressed that the 
negative value of the derivative dSciXiQ 2 ) / das is due to the small x contribution to the 
Gottfried integral. 

Figure 2 shows the dependence of Sg(Q 2 ) on asj^ at x = 10~ 3 in (7). The in- 
dependence could be parametrized by the parabola: 

S G (Q 2 ) = S ■ (1 + c x ■ [ptsM + c 2 ■ K/vr) 2 ). 



The values of S , ,Ci,c 2 are presented in Table 1. They could be considered as as 
corrections to the Gottfried sum rule. One can see from Figure 2 that the derivative 



(1Sg{Q 2 ) / dots is negative at as/ft < 0.05 for all parametrizations F 2 , and the HI curve 
crosses the NMC one in the range Q 2 = 5 — 10 (GeV/c) 2 . 



Table 1. The coefficients of the Sg(Q 2 ) Gottfried integral expansion in a^/vr. 





So 


Cl 


Cl 


NMC92 


0.246 


-7.06 


53.2 


NMC95 [| 


0.210 


-4.87 


44.2 


HI [11 1 


0.271 


-6.04 


23.7 


average value 


0.242± 0.21 


-6.00± 0.74 


40.4± 11.1 



Despite different kinematical regions of NMC and HI experiments used for fits of 
experimental data one can see from Table 1 the reasonable quantitative and good quali- 
tative agreement between the corresponding coefficients So, c±, c 2 for parametrizations of 
the proton SF under consideration. Note that the F 2 parametrization of HI is combined 
with the data from NMC and BCDMS experiments and a smooth transition between dif- 
ferent Ff data is obtained. The kinematic range covers almost four orders of magnitude 
in x and Q 2 . 

The average values for the coefficients (i = 5 , ,c 1 ,c 2 ) are calculated by < r >= 
{J2i^i)/ n and < Ar >= (Yh\ < r > ~ r i\)/ n an d presented in Table 1. The errors 
obtained could be considered as a crude estimation of systematic uncertainties. 

It should be noted that the parametrization of the ZEUS data |12| for F 2 provides 



qualitatively the same results: So = 0.383, c\ = —12.9 and c 2 = 76.2. 

We would like to emphasize that the values of coefficients obtained from the phe- 
nomenological analysis of experimental data are essentially different from the theoretical 
QCD predictions for c\ and c 2 . The coefficient c\ is found to be negative, in contrast to 
c\. Both ci, c 2 are many times larger in the absolute value than c±, c 2 . 



7 Conclusions 

The analysis of Q 2 -dependence of the Gottfried sum rule based on experimental data 
on the proton F 2 and deuteron F 2 D structure functions in the framework of the covariant 
approach in the light-cone variables and relativistic deuteron model was performed: 

• the procedure to extract the neutron SF F£(x, Q 2 ) is described and used to analyze 
the parametrizations of NMC, BCDMS, SLAC, HI, ZEUS data on F 2 D and Ff SF's 

• the increase of the Sg(Q 2 ) Gottfried sum rule for small o>s{Q 2 ) is a general fea- 
ture for the NMC92, NMC95, HI parametrizations (see Fig. 2). This behaviour is 
connected with the negative value of the first order as correction to the GSR 



• it is shown that the results obtained for c\ and C2 are in disagreement with the 
calculation made in the framework of QCD assuming the flavour symmetry of sea 
quarks (u = d) 

• the measurements of the deuteron structure function F 2 D at HERA, extraction the 
neutron SF i 7 ^ 1 and verification of the Q 2 - dependence of Sq(x, Q 2 ) at low x are 
necessary to determine the a^-corrections to GSR more reliably. 

Acknowledgement 



The authors would like to thank Yu.A.Panebratsev for his support of our work. This 
work was partially supported by Grants of the Russian Foundation for Fundamental 
Research under No. 95-02-05061 and No. 96-02-18897. One of us (M.V.T.) are very 
grateful to R.Klanner and M.Klein for support to take part in the Workshop "Future at 
HERA". 



References 

[1] K. Gottfried, Phys. Rev. Lett. 18 (1967) 1174. 

[2] J. Bjorken, Phys. Rev. 148 (1966) 1467. 

[3] J. Ellis, R. Jaffe, Phys. Rev. D9 (1974) 1444. 

[4] D. J. Gross and C. H. Llewellyn Smith, Nucl. Phys. B14 (1969) 337. 

[5] S. L. Adler, Phys. Rev. 143 (1966) 1144. 

[6] P. Amaudruz et al., Phys. Rev. Lett. 66 (1991) 2712; 
M.Arneodo et al, Phys.Rev. D50 (1994) Rl. 

[7] P. Amaudruz et al., Preprint CERN-PPE/92-124, 1992. 

[8] M.Arneodo et al, Phys. Lett. B364 (1995) 107; 
Preprint CERN-PPE/95-138, 1995. 

[9] M.R. Adams et al, Phys. Lett. B309 (1993) 477. 

[10] M.Diiren, DESY HERMES-95-92, 1995. 

[11] T.Ahmed et al, Nucl. Phys. B439 (1995) 471; 
A.Aid et al., DESY 96-039, 1996. 

[12] M.Derrick et al., Z.Phys. C65 (1995) 379. 

[13] M.A.Braun, M.V.Tokarev, Particles and Nuclei 22 (1991) 1237. 
[14] M.A.Braun, M.V.Tokarev, Phys. Lett. B320 (1994) 381. 



[15] D.A. Ross, C.T. Sachrajda, Nucl. Phys. B149 (1979) 497; 

F.J.Yndurain, Quantum Chromodynamics (An Introduction to the Theory of 
Quarks and Gluons).- Berlin, Springer- Verlag (1983), 117. 



[16] A.L.Kataev et al., Preprint INP-0919/96, Moscow, 1996 ||hep-ph/9605367 
Submitted to Phys. Lett. B. 

[17] L.Whitlow et al., Phys. Lett. B282 (1992) 475. 

[18] A.C. Benvenuti et al, Phys. Lett. B223 (1989) 485; 
A.C. Benvenuti et al, Phys.Lett. B237 (1990) 592. 

[19] B.Badelek, J.Kwiecinhski, Nucl. Phys. B370 (1992) 278. 

[20] A.V. Sidorov, M.V. Tokarev, Phys. Lett. B358 (1995) 353. 



a) 



b) 



c) 



d) 



e) f) 

Figure 1. The Gottfried integral Sq(x, Q 2 ) and the derivative dScix^Q 2 ) / das as a 
function of x and Q 2 . 



Figure 2. The Gottfried integral Sg{,Q 2 ) as a function ofa:,sat:r = 10 3 . The 
lines present the paraboloic fit of Sa(, Q 2 ) for different parameterisations F%. 
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